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ABSTRACT 

NGC 4753 is a bright (My ~ 22.3) lenticular galaxy. It is a very interesting target to test 
different theories of formation of lenticular galaxies, due to its low-density environment and 
complex structure. We perform the first comprehensive study of NGC4753 globular cluster 
system (GCS), using Gemini/GMOS and CTIO/MosaicII images. Our results indicate a rather 
poor GCS of Ri 1000 members. Its azimuthal distribution follows the shape of the galaxy 
bulge. The GC colour distribution is peculiar, presenting an intermediate subpopulation in 
addition to blue and red ones. This intermediate subgroup can be explained by a single stellar 
population with an age of 1.53 Gyr and 0.5 — 1 Zq. The GC specific frequency Sn = 1-3 ± 
0.15 is surprisingly low for a galaxy of its class. The GC luminosity function (GCLF) is 
also peculiar, with an excess of bright GCs compared to the expected gaussian distribution. 
The underlying galaxy body has significant substructure, with remnants of spiral arms, dust 
filaments, and isophote twisting. This, and the fact that NGC4753 hosted two type la SNe, 
support the possibility that the intermediate GC subpopulation may have originated during a 
recent merger, 13 Gyr ago. 

Key words: galaxies: star clusters: individual: NGC 4753 - galaxies: elliptical and lenticular, 
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1 INTRODUCTION 

Despite that the fraction of lenticular ga laxies (SOs) seems to in - 
crease with the number of neighbours l lWilman & Erwinll2012h . 
the existence of SOs in low density environments poses ques¬ 
tions about their origin. The physical processes involved in the 
evolution of lenticula r galaxies (SOs) have b e en largely discussed 
in the literature (e.g. ISnitzer & Baadelll95ll : I Gunn & God 1 197 2l: 


__^en Bergh|[^09hh . Ram pr essure stripping (e. g. Moran et all 
20061) and galaxy harrassment jMoore et alJll99^ are currently 


considered as plausible explanations for the origin of SOs in 
high density environments. However, these processes are not ex¬ 
pected to play a main role in low density environments. Alterna¬ 
tive explanations are disk gas removal by an acti ve galactic nu¬ 
cleus, or an inefficien t inflow of external gas (e.g. iLi et alJ[200^: 
Ivan den B ergh||2009a|. a nd references therein). How ever, cases like 
NGC 404 lldel Rio e(alj[200i: iThilker et al.ll2010l) and NGC 4460 
jMoiseev et al .120101) seem to present a steady, long term gas accre¬ 
tion, with current star formation, implying a complex evolutionary 
history. 


* E-mail: jpcaso@fcaglp.unlp.edu.ar (JPC), lbassino@fcaglp.unlp.edu.ar 
(LPB), matiasgomez@unab.cl (MG) 


NGC 4753 is an early-type galaxy (ETG) with prominent and 
twisted dust lan es. It is considered as the brightes t galaxy of one of 
the MK-groups jKarachentsev & Nasonovall2013l) that populate the 
Virgo S outhern Extens ion, which is assigned 23 members. Alter¬ 
natively, ) ll993l) has included NGC 4753 in the one dozen- 

members’ group EGG 315 (or the ’NGC 4643 group’ according to 
NED, NASA/IPAC Extragalac tic Database). O riginally classified 
as a peculiar lenticular galaxy l lSandagelll96ll) . it is now assigned 
an irregular de Vaucouleurs type 10 (e.g. NED). This latter clas¬ 
sification points to a merger or a galaxy distorted by gravitational 
interaction with a close neighbor, though the existence of a neigh¬ 
bor is not evident. The apparently high type la Supernovae (SNIa) 
rate exhibi ted by this galaxy (N G C 4753 hosted SN 19651 and 
SN 1983G.ICiatti & Barbo^ll97ll and iKosai et alJl983l:lButa et all 
1 19851 : iMueller & Hoeflichlll994 respectively) may be an indica¬ 
tion of an intermediate-age population, a stellar population that 
may have formed in a merger a few Gyrs ago. However, this 
latter statement seems c ontradictory to the study performed by 
iNavasardvan et all (1200 ll) . who concluded that excepting strongly 
interacting systems, the SN production in isolated galaxies, galaxy 
pairs and groups does not depend on the environment. Or has 
NGC 4753 possibly suffered such a strong event?, we will come 
back to this issue in the Discussion. 
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Regarding the dust content, ISteiman-Cameron et alj l ll992t) 
proposed a disk model that is fitted to the dust distribution, so that 
the intrincated lanes can be understood as the result of a disk that 
has been twisted deeply by differential precession. In this way, the 
original S O turns into a peculiar galaxy after an accretion event. 
Moreover, IPewangan et alj d 19991) have estimated photometrically 
the mass of cold dust with optical and far-infrared (FIR) data, and 
showed that it is about a factor 10 higher than t hat estimated by 
mass-l oss from only red giant stars. More recentlv. fFinkelman et al.l 
(l2012h studied the relation between dust and ionized gas in ETGs 
with dust lanes, including NGC 4753 in their sample. In order to ex¬ 
plain the origin of the interstellar medium (ISM) detected in these 
ETGs, in addition to the internally-produced ISM, they stressed 
on the plausible capture of a gas-rich spiral or dwarf, or the col¬ 
lision between two gas-rich dwarfs with similar masses. Thus, we 
are again led to an accretion or merger event. 

With the aim of disentangling the evolutionary story of 
NGC 4753 and getting a new and precise distance determination, 
we present a deep photometric study of its globular cluster (GC) 
system. It is well known that old globular cluster systems (GCS), 
i.e. with ages greater than lOGyr, carry a historical record of the 
star formation eve nts and, more generally, of the evolution of their 
host g alaxies (e.g. lBrodie & Stradeill200a : lToninil|2013l : [Kniiisseiil 
|2014 and references therein). The formation of massive star clus¬ 
ters happens during major star-bursts. Thus, the GC colour distribu¬ 
tion can provide evidence of star-burst episodes or merger events, 
which can be as sociated with the pre s ence of an interme diate-age 
popula tion (e.g. iR ichtler e t al.ll201 j ; ICaso et alj 120139) . Eor in¬ 
stance, ISahnas^et^ 12015 ) have indicated that low v alues of the 
specific frequency Sn toarris & van den BerghifIgsTI) or Tjv, be¬ 
ing Tjv the number of GCs normal ized by the stellar mass of the 
host galaxy jZepf & Ashmaniri993h . are due to the paucity of red 
(bona fide high metal content) GCs. In the current scheme of GCs 
origin, blue (bona fide low metal content) GCs are formed at high- 
redshift, during the massive star formation episodes that follow the 
merging of the first building blocks of proto-galaxies. Subsequent 
merger events are responsible for the formati on of red GCs, ac¬ 
companied by another bu rst of few blue ones l Muratov & GnedinI 
120101 : IlI & Gnedii]|2014h . Alternatively. iToninil ( 1201 3I) ' presents a 
hierarchical clustering model where the metal-rich subpopulation 
is composed of globular clusters formed in the galaxy main pro¬ 
genitor around redshift z ~ 2, and the metal-poor subpopulation 
is composed of clusters accreted from satellites, and formed at red- 
shifts 2: « 3 — 4. In this context, the study of GCSs will help to 
understand the origin and evolution of ETGs in low density envi¬ 
ronments. 

Distance determinations for NGC 4753 range between 20 


and 2 4 Mpc, including sup ernova light-curve analysisJParodie^^ 


^_ up i _ 

200c; iReindl et'H] 1 20051) . Tully-Eisher r elation j Theureau et al.. 
2007), and surface brightness fluctuations jToni2.et^ TOOjJ). 


will adopt a mean distance of « 23.6 Mpc jTonrv et alj200l[) in the 
rest of this paper, until we can compare it with the result obtained 
from the GC luminosity function (GCLF). There are no previous 
distance measurements using the turn-over of the G CLF. The infor¬ 
mation on NGC4753’s GCs is rather non-existent. IZaritskv et al.l 
( l2015l) include NGC 4753 in their sample of 97 ETGs for which 
the number of GCs as well as radial profiles are estimated using 
Spitzer S4G data. However, as the authors indicate, these are rather 
crude estimates and focused on the central region. 


This paper is organized as follows. The Observations are pre¬ 
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Figure 1. Positions of the two GMOS-S fields from our programme, su¬ 
perimposed to an R image from the Palomar Observatory Sky Survey. The 
image size is 15' X 7.5'. North is up, East to the left. 


sented in Section 2, the Results are shown in Section 3, and the Sec¬ 
tion 4 is devoted to the Discussion. 


2 OBSERVATIONS AND REDUCTION 
2.1 Observational data 

Our dataset consists of GMOS-Gemini-South images in g', r’, 
and i' filters, observed during March/April 2014 (programme 
GS2014A-Q-67, PL J.P. Caso). These images correspond to two 
fields, one of them centred on NGC 4753 (hereafter, N4753F), and 
the other was located to the West (hereafter, WestF). Both fields 
were slightly overimposed, in order to determine possible differ¬ 
ences in the zero-magnitude (see FigO. The exposure times were 
4 X 570s in g', 4 x 270s in r', and 4 x 360s in i'. Each set of expo¬ 
sures were slightly dithered, in order to fill the gaps in the GMOS 
field, and to efficiently remove cosmic rays and bad pixels. The re¬ 
sulting seeing for the three fi lters was 0.6" — 0.8". The a ssumed 
absorptions come from NED jSchlaflv & Finkbeineill20ll]) . 

The foreground contamination by Galactic stars was esti¬ 
mated from the observations of the William Herschel Deep Eield 
(WHDF). This was necessary bevause the GCS extended over 
the two science fields (see Section lT^ . The reduced observations 
of this field, part of the GMOS-N system verification observa¬ 
tions (programme GN-201 lB-SV-104), were downloaded from the 
Gemini Science Archive. This field had been observed in the same 
filters than our programme, with similar exposure times. Accord¬ 
ing to its Galactic Latitude (b = —61.7°), the foreground contam¬ 
ination would be similar to that of the N4753E (h = 61.7°). The 
reduction was perfomed in the usual way, using the tasks from the 
GEMINI-GMOS package, within IRAF. It was necessary to sub¬ 
tract the fringe pattern in the i' images from our programme. The 
bias, flatfield, and i ’ blank sky images were obtained from the Gem¬ 
ini Science Archive. 

This dataset is further enhanced by using images in Kron- 
Cousins R filter, obtained with the MOSAIC II camera mounted 
at the 4-m Blanco telescope at the Cerro Tololo Inter-American 
Observatory (CTIO, Chile). The observations were carried out on 
March 25-26, 2004. The observation was split in four exposures of 
720 s each. During these nights, one of the MOSAIC II CCDs was 
not working. To overcome this, individual images were taken with 
an appropiate dither to fill in the missing chip. Thus, the final image 
presents a tield-of-view of 36 x 45 arcmin^, i.e., « 247 x 309 kpc^ 
at the assumed distance to NGC 4753. 
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2.2 Photometry and point source selection 

First, the extended galaxy light was substracted from the CMOS 
obser vations of the N4753F and WestF, applying a ring median 
filter ( ISeckeiiri995h with an inner radius of Sarcsec (35 px), and 
an outer radius of barcsec (41 px). This procedure facilitates the 
point source detection, and is commonl y used in GCs studies (e.g . 

IPirsch et al.l2003l : lBassino et al.l2006al la: l&so et alj2013aLl201^ . 

This procedure was not necessary in the case of the WFID F. 

The software SEXTRACTOR jBertin & Amoutsll996h was ap¬ 
plied to the i' image (because it presented the higher signal-to- 
noise ratio), in order to obtain an initial point source catalogue. 
The effective radii of clasi cal GCs is usually smaller than 10 pc 
(e.g. lBriins & KroupJ20121) . Considering that, at the assumed dis¬ 
tance, 1" corresponds to « 110 pc, the NGC4753 GCs are seen as 
point sources on our images. Then, we used the SEXTRACTOR pa¬ 
rameter CLASS_STAR to eliminate the extended sources from our 
catalogue. 

The phot ometry was performed with the DAOPHOT package 
dStetsonl 19871) within IRAF. A second-order variable PSE was gen¬ 
erated for each filter, employing a sample of bright stars, well dis¬ 
tributed over the field. The final point source selection was made 
with the and sharpness parameters, from the ALLSTAR task. 
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Figure 2. Completeness curves for the three fields (our science fields plus 
the WHDF), as function of Iq magnitude. The bin width is 0.25 mag. 


within the uncertainties. Considering that, it was possible to calcu¬ 
late a single set of transformation equations, which resulted 


C = c-0.059±0.004-0.418±0.010xAc-F0.111±0.005x(C'-ri) 

Ti = r+0.628±0.005-0.14±0.002xAfl+0.019±0.002x(C'-Ti) 

where r and c are the instrumental magnitudes, and Xr and Xc 
are the airmasses in R and C filters, respectively. 


2.3 Photometric calibration 

During o ur programme, a wei l popuiated standard star field, from 
the list of ISmith & et alj | |2002h was observed during the same night 
than the N4753F. We obtained the standard stars aperture photom¬ 
etry for several aperture radii with task PHOT. We considered these 
measurements to derive the aperture correction by using a growth 
curve. After this correction, we fit transformation equations of the 
form: 

mstd = ZP -F rriinst - Kcp x (X — 1) -F CT x (mi — m 2 )atd 

where rristd and rriinst are the standard and instrumental magni¬ 
tudes, respectively, and ZP are the photometric zero points. Kcp 
is the mean atmospheric extintion at Cerro Pachon, obtained from 
the Gemini Observatory Web Pag^ and X the airmass. CT is the 
coefficient of the color term, and (mi — m 2 )atd its associated cal¬ 
ibrated colour, {g' — r') for g' and r', and (r' — i') for i'. The 
fitted zero points and colour terms were ZP® = 28.185 ± 0.013, 
ZP’~' = 28.314±0.009 and ZP*' = 27.915±0.010, and CTg/ = 
0.106±0.018, CTr> = 0.070±0.012 andCTi/ = 0.039±0.023, 
respectivel y. The calibration equations for the WHDF were ob¬ 
tained from iFaifer et al.l l l201lh . 

Afterwards, we applied the galactic extinction corrections 
from ISchlaflv & Finkbeinej ( 1201 ll) to the calibrated magnitudes. 
Finally, from the objects in common between the N4753F and 
WestF we found small zero-point differences, Agi = —0.05, 
— 0.04 and A^/ = 0.02. These offsets were applied to 
the WestF final catalogue, in order to refer its photometry to the 
N4753F one. 

In the case of the MOSAIC II observations, 4 to 5 fields were 
observed in Kron-Cousins R and Washi ngton C filter s , conta ining 
about 10 standard stars from the list of iGeisler et al.l l ll996l) . The 
observations spanned a large range in airmass, typically from 1.0 
to 1.9. The coefficients fitted for each night were indistinguishable 

^ http://www.gemini.edu/sciops/instruments/gmos/calibration 


2.4 Completeness analysis 

In order to estimate the photometric completeness for the three 
fields (our science fields plus the WHDF), we added 250 artificial 
stars in the images of the three filters. These stars were equally dis¬ 
tributed over the entire GMOS-S fields, they presented colours in 
the expected ranges for GCs, and i > 20.75. This process was re¬ 
peated 40 times, achieving a sample of 10000 artificial stars. Their 
photometry was carried out in the same way as the science fields. 
The completeness curve, after the definite point source catalogue 
was generated, is shown for both fields in Figure[2 The complete¬ 
ness functions are very similar, achieving the 80% at i! « 24.1. 
This value was considered as the faint magnitude limit in the fol¬ 
lowing analysis. 

The comparison between the initial brightness of artificial 
stars and their subsequent photometry proved that the application 
of the median filter does not modify the point source photometry. 

3 RESULTS 

3.1 GC candidates selection from colours 

Figure[3]shows {g' — T)o vs. {r' — i')o and (g' — r')o vs. (g' — i')o 
colour-colour diagrams (CCD) for point sources in both fields. Old 
GCs present narrow colour ranges, which can b e helpful to distin¬ 
guish them from cont amination in the CCDs (e.g. lFaifer et alj201 ll : 
lEscudero et alj 1201^ . It is particularly evident in the N4753F, 
where it is expected to find the majority of the GCs belonging to 
NGC4753 system. We selected as GC candidates the point sources 
ranging 0.4 < (g' — i')o < 1.4, 0.3 < {g' — r')o < 1 and 
0 < (r' — i')o < 0.5. These values are simi lar to those chosen by 
lFaiferetal.l ( l201lh and lEscudero et all ( l2015h . The number of point 
sources brighter than i' — 24.1 (our 80% completeness limit, see 
above) that fulfill the colour criteria are 437 in the N4753F, 117 in 
the WestF, and 32 in the WHDF. 

In Figure|4]we show the ig vs. {g' — T)o colour-magnitude 
diagrams for the three fields. The black dots represent the point 
sources, after the x^ and sharpness selection. The blueviolet filled 
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Figure 3. Colour-colour diagrams for both science fields. Solid lines indi¬ 
cate the colour ranges of GC candidates (see Section 3.1). 
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Figure 4. Colour-magnitude diagrams for the three fields. Black dots repre¬ 
sent point sources, while blueviolet filled circles highlight those objects that 
fulfill the colours and magnitude criteria applied to select GC candidates. 


circles highlight those objects that fulfill the colours and magni¬ 
tude criteria indicated before. The CMD in the middle panel, corre¬ 
sponding to N4753F, shows that GC candidates span a rather large 
range, with a large fraction of blue GCs. Comparing the CMDs 
from the WestF and WFIDF, it can be seen that some GC can¬ 
didates are still present in the former one. If we assume that the 
brightness limit between ultra-compact dwarf s (UCDs) and ‘clas- 
sical’ GCs is My ~ —10.5 — —11 fe.g. iMieske et al.l l2006l : 
lHilkeill2009l: ICaso et al.ll2013aLl2014h . at NGC4753 distance im¬ 
plies Iq « 20.9 — 20.4. Hence, NGC 4753 does not seem to present 
any candidate to UCD. 


3.2 Colour distribution 

The Figure|5] shows the colour distribution for the GC candi¬ 
dates with ig < 24.1, applying a binwidth of 0.04 mag. The 
colour distribution seems to deviates from the ty pical bimodal pat¬ 
tern, frequently se en in early-type galaxies (e.g. lForte et al. I |2014 
ISrodie et al.ll2014l and references there in). We fit means and dis¬ 
persions for bimodal and trimodal distributions, us ing the Gaussian 
Mixture Method IGMM. lMuratov & Gnedinl201(]h . The results are 
listed in TableQ] where we also include the fractions of each sub¬ 
population / as well as the DD parameter, that indicates whether 
an specific multimodal distribution provides a realistic fit (where 
the null hypothesis is accepted for DD> 2). For the three com¬ 
ponents distribution, the single Gaussians are plotted with dashed 
lines. Doted lines represent the Gaussians for the bimodal distribu¬ 
tion. 

If we consider all GC candidates together (first and second 
rows in Table[T}, we can see that for the trimodal distribution, 
the bluer and redder Gaussian means resemble those of blue and 
red GC (metal-poor and metal-rich) subp opulations, respectively 
(e.g. lFaifer et aljEoill : I Salinas et al.ll2015l) . However, there is evi¬ 
dence of a third population, presenting intermediate colours. This 
may reflect a complex evolutionary history, si milar to other early- 
type galaxies in low- density environments (e.g. ICaso et aFlEoBbl: 
lEscudero et alj2015h . The intermediate GC subpopulation presents 
mean(y' —i')o « 0.9 in all cases (TablellJ. We compare this colour 
with the theoretical models of single stellar populations (SSPs) by 
ISressan et al.l ( l2012h . using their web-bas ed toolEl If we consider 
SSPs with metallicities 0.5 Zq and Zq, a lChabri^ ( 1200 ih lognor¬ 
mal IMF, and absorption Ay = 0.09 (NED), the mean colour 
for the intermediate subpopulation corresponds to a population of 
1.5 — 3Gyr. This age is consistent with the existence of SNIa in 
the galaxy population. The derived ages shoul be interpreted as raw 
estimations. The selection of these particular metallicities, as well 
as the IMF, might introduce systematic errors. 

Another interesting issue is that if we consider the bimodal 
distribution, the mean value, dispersion and fraction of blue GCs 
are in agreement, within uncertainties, with the equivalent values 
of blue GCs from the trimodal distribution. However, the fraction 
of red GCs in the bimodal case, seems to include approximately 
the intermediate plus the red subpopulations of the trimodal case. 
That is, the blue subpopulation remains almost unchanged. This 
might be the reason for the overestimation of the red population in 
Figure|5] 

In order to look for radial gradients, we separate the sample 
in two radial regimes. We assume an elliptical projected spatial 
distribution, with PA and ellipticity equal to 91° and 0.45 (see 
Section[33). The adopted radial galactocentric limit is taken as 
Rg = 2.5 arcmin, which results in two samples with similar num¬ 
ber of GC candidates. The results for both radial ranges are de¬ 
picted in the four following rows of Table[T] In both cases, the DD 
values from GMM indicate that a trimodal profile could represent 
the colour distribution. The means and dispersions are similar, but 
the intermediate subpopulation dominates the distribution when we 
move towards larger galactocentric distances, at the expense of the 
blue subpopulation. 

The comparison between faint and bright GCs also reveals 
different behaviours in the colour distribution. The selection of 
ig = 23 as the adopted magnitude limit to separate bright and faint 
GCs is based on the GCs CMD (Fig.|4l(, and luminosity function 
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Figure 5. Upper panel: Colour distribution for the GC candidates. The 
dotted lines indicate the Gaussian fits to the bona fide metal-poor popula¬ 
tion and metal-rich population assuming a bimodal distribution. For three 
components, the Gaussian fits are shown with dashed lines. The solid lines 
represent the sum of Gaussians for each scenario. Lower panel: Colour 
distribution for GCs candidates in the magnitude range 24.1 > Iq > 23 
(dashed line histogram) and Iq < 23 (solid line histogram). 


(GCLF) (see Section lT!4l Despite of the similarity in means and 
dispersions, the weight of the different populations in the colour 
distribution differs significantly. The faint GCs seem to present a 
trimodal distribution, as well as the entire sample. The presence 
of the red population in the bright sample seems marginal (lower 
panel of Fig[3- 

Thus, both bimodal and trimodal distributions seem to repre¬ 
sent the colour distribution with different degree of accuracy. Flow- 
ever, the dependences of the distribution with galactocentric radius 
and brightness point to a GCS built with three different compo¬ 
nents. 


3.3 Azimuthal and radial distributions 

The upper panel of Figure[^shows the azimuthal distribution for all 
the GC candidates at less than 2.5 arcmin from the galaxy centre. 
This limit is chosen to achieve an homogeneous azimuthal cover¬ 
age, according to the field-of-view of the CMOS field. The posi¬ 
tion angle (PA) of the GCs projected distribution is measured from 
North through East, and the bin width is 20°. Considering that an 
elliptical GCS causes a sinusoidal distribution in this diagram, we 
fit a cosine function to the histogram. The fitted PA is 94±6°, which 
is in agreement with the galaxy light PA (see Section iTSt . The am¬ 
plitude and the symmetry axis offset are 4.7 ± 0.9 and 8.9 ± 0.7, 
respectively. 

The lower panel of Figure|6]exhibits the logarithm of the back¬ 
ground corrected radial distribution for the GC candidates. The 



Figure 6. Upper panel: azimuthal distribution for GC candidates at less 
than 2.5 arcmin from the galaxy centre. The GCs position angle (PA) is 
measured from North through East, and the bin width is 18°. Lower panel: 
radial distribution for GC candidates. The solid line represents the power- 
law fitted by least-squares. The dashed line coresponds to the background 
level, while the dotted one indicates the 30% of the background level, which 
is the used to define the GCS extension. 


population of the GCS is not large enough to discriminate between 
different subpopulations. 

Due to the elongated spatial distribution of the GCS, we con- 
sidere concetric ellipses instead of circles. The ellipses PA is the 
previously fitted one, while the ellipticity was assumed as « 0.45, 
obtained from the galaxy diffuse light (see Section[33J. The hor¬ 
izontal dashed line indicates the background level, while the dot¬ 
ted one is the 30% of it. We propose this limit to define the GCS 
extension, which has been previously used in similar studies (e.g. 
ISassino et al.l200^ : ICaso et alj2013^ . The background corrected 
distribution can be fit by a power-law, with slope —3.2 ± 0.3. Ac¬ 
cording to this, the total extension of the GCS is « 8 arcmin (i.e., 
« 50 kpc). The inner bin of the radial distribution seems to under¬ 
estimate the GCs density. This is probably due to a higher imcom- 
pleteness, caused by the high surface brightness of the galaxy, plus 
its striking dust structure. 

The projected spatial distribution for GC candidates for two 
different magnitude bins is shown in Figure^ The colour gradient 
in both panels spans the range 0.4 < (g' — i')o < 1.4. North is 
up, and East to the left. It seems that GC candidates brighter than 
i'o = 23, mainly blue candidates, present inhomogeneities in their 
spatial distribution. In particular, there is a paucity of GCs to the 
North of NGC4753. The faint GC candidates present a much ho- 
mogeneus spatial distribution, although they follow the major axis 
of the galaxy (approximately in the East-West direction). It can be 
seen GCs from the three subpopulations at large projected distances 
from the the galaxy. This deviates from the usual scenario, where 
the b lue GCs are mor e spread, following t he X-ray gas emissio n 
(e.g. iForte et all l2005l : iForbes et"^ 1201A lEscudero et al.l 120151) . 
and the red ones are more concentrated, following the galaxy light 
(e.e. lSassino et alj2006ah . 
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Table 1. Parameters of the GMM fitting to the colour distribution, considering bimodal and trimodal cases, fii, <Ji and fi coiTespond to the means, dispersions 
and fractions for each Gaussian component. The parameter DD indicates whether an specific multimodal distribution provides a realistic fit (where the null 
hypothesis is accepted for DD> 2). 





h 

M2 

cr2 

h 

M3 

0-3 

h 

DD 

All GC candidates 

Trimodal 

0.76 ±0.05 

0.11 ±0.02 

0.63 ±0.13 

0.91 ±0.04 

0.03 ± 0.02 

0.11 ±0.05 

1.05 ±0.02 

0.07 ±0.01 

0.25 ±0.06 

2.04 

Bimodal 

0.80 ±0.04 

0.13 ±0.02 

0.81 ±0.15 

1.06 ±0.05 

0.07 ±0.02 

0.19 ±0.15 




2.50 

Rg < 2.5' 

Trimodal 

0.76 ±0.05 

0.11 ±0.02 

0.76 ±0.07 

0.91 ±0.03 

0.03 ±0.01 

0.13 ±0.06 

1.04 ±0.02 

0.05 ±0.01 

0.11 ±0.02 

2.15 

Bimodal 

0.80 ±0.03 

0.12 ±0.01 

0.95 ±0.04 

1.06 ±0.05 

0.04 ± 0.02 

0.05 ±0.04 




2.85 

Rg > 2.5' 

Trimodal 

0.69 ±0.05 

0.09 ± 0.03 

0.33 ±0.10 

0.89 ±0.04 

0.06 ± 0.03 

0.31 ±0.12 

1.06 ±0.03 

0.07 ±0.01 

0.36 ±0.09 

2.60 

Bimodal 

0.68 ±0.08 

0.08 ± 0.03 

0.41 ±0.17 

0.96 ±0.05 

0.10 ±0.03 

0.59 ±0.17 




2.55 

ig < 23 

Trimodal 

0.74 ± 0.04 

0.07 ±0.02 

0.54 ±0.08 

0.92 ±0.04 

0.05 ±0.02 

0.40 ± 0.06 

1.14 ±0.05 

0.03 ±0.01 

0.06 ±0.01 

2.92 

Bimodal 

0.69 ±0.06 

0.15 ±0.04 

0.44 ±0.14 

0.87 ±0.13 

0.06 ±0.01 

0.56 ±0.14 




2.10 

23<i'g< 24.1 

Trimodal 

0.75 ±0.03 

0.11 ±0.02 

0.57 ±0.10 

0.91 ±0.02 

0.04 ±0.01 

0.14 ±0.07 

1.06 ±0.01 

0.06 ±0.01 

0.29 ±0.04 

2.20 

Bimodal 

0.80 ±0.03 

0.13 ±0.02 

0.77 ±0.09 

1.06 ±0.03 

0.06 ± 0.02 

0.23 ±0.09 




2.62 
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Figure 7. Projected spatial distribution for GC candidates, split in two mag¬ 
nitude regimes. The centre of the galaxy is highlated with a cross. North is 
up, and East to the left. The field of view is 9.1 X 5.5 arcmin^. The colour 
gradient spans the range 0.4 < {g' — i')o < 1.4. 


3.4 Luminosity function and total nnmber of GCs 

In order to facilitate the calculus of the total number of GCs, we ob¬ 
tain the GCLF for those GCs located within two concentric ellipses 
with the same PA and ellipticity, equal to 91° and 0.45, respectively 
(see Section[33J. major semiaxes 1.25 and 3.6 arcmin. Figure[8] 
shows the raw GCLF with open circles, and the background and 
completeness corrected one with blueviolet filled circles. The error- 


bars assume Poisson uncertainties for science and background mea¬ 
surements, and the binwidth is 0.2. The vertical lines indicate the 
luminosity range that has not been considered in the GCLF fitting, 
due the declining completeness (we selected the limit Iq = 24.1 
in Section 2.4). It can be noticed a change in the GCLF behaviour 
around iq « 23, showing an excess of bright GCs in comparison 
with faint ones. We fitted a Guassian profile to the entire GCLF, 
which is represented by the dotted curve. The resulting turn-over 
magnitude (TOM) and dispersion are ig tom ~ 24.2 ± 0.2 ± 0.1 
and a = 0.8 ± 0.12 ± 0.1, where the second uncertainty comes 
from the binning considered here. If we fit just the GCLF bins with 
i' > 23, we obtain a Gaussian profile slightly different, which is 
plotted with a dashed curve. Comparing both profiles, it is clear 
that the excess of bright GCs results in a higher dispersion for 
the fitted Gaussian, and consequently the TOM moves slightly to¬ 
wards fainter values. The TOM and dispersion in this latter case 
are Iq^tom = 24 ± 0.08 ± 0.1 and cr = 0.59 ± 0.07 ± 0.1, re¬ 
spectively. In the following, we will assume that these parameters 
represent the GCLF for NGC 4753. 


Usually, old GC populations present a Gaussi an GCLF, with 
a TOM in the U-band of Mtom ~ —7.4 (e.g. lRichtleill200^ : 
llordan et al.l2007l) . The GC candidates fainter than ig = 23 present 
a mean colour {g' — i') « 0.9 ± 0.03, and then g'n t o m ~ 24.9 ± 
0.12. Applying Equation 2 from iFaifer et alj j201lh . the TOM in 
the V filter is Vq^tom ~ 24.55 ± 0.14, which implies that the dis¬ 
tance moduli and metric distance result (m — M) ~ 31.95 ± 0.14 
and « 24.5 ± 1.5 Mpc, respectively. This distance is in agree ment 
with previous mea surements from S BF jTonry et al.l 200 ih and 
SNIa light-curves jParodi et aP 120001: iReiiidr et a D l2005h . GCLF 
dispersion seem to be lower than expected. 


The numerical integration of the Gaussian profile results in 
445 ± 35 GCs. We have to consider that the GCLF was obtained 
from the GC candidates in a specific radial regime, that does not 
span the entire GCS extension. From the density profile derived 
in Section[33] we calculate the ratio of GCs in the radial regime 
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Figure 8. Luminosity function (GCLF) for GCs in the radial regime 
1.25 — 3.6arcmin (assuming an elliptical geometry). Open circles repre¬ 
sent the raw GCLF, while blueviolet filled circles indicate the background 
and completeness corrected one. The errorbars assume Poisson uncertain¬ 
ties for science and background measurements, and the binwidth is 0.2. The 
vertical lines indicate the luminosity range that has not beenconsidered in 
the GCLF fitting. 
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1.25 — Sarcmin to GCs in the radial regime 1.25 — 3.6arcmin, 
« 1.26 ± 0.07. Due to the incompleteness, we cannot derive the 
surface density of GCs in the inner 1.25arcmin. A common cri¬ 
terion in GC studies is to adopt for this region the same surface 
density as for the first bin (e.g. ISassino et alj|200^ : ICaso et ^ 
l2013bh . and correct it for completeness. This can be done, assum¬ 
ing that the radi al distribution of GCs flattens in the inner region 
of a galaxy (e.g. lElson et IDl 19981) . Hence, we obtain 1030 ± 120 
GCs. We have not considered yet the excess of bright GCs. We can 
do that subtracting to each GCLF bin the numerical integration of 
the Gaussian profile fitted to the GCs fainter than i' = 23. The sum 
of this excess result in « 40 GCs. Then, the number of members 
of the NGC4753 GCS is 1070 ± 120. Assuming the tota l V mag¬ 
nitud e obtained by the Camegie-Irvine Galaxy Survey ( IHo et al.l 
lioTII) , Vtot = 9.66 ± 0.11, and the distance moduli previously 
derived, m — M = 31.95 ± 0.14, the absolute magnitude of the 
galaxy is My = —22.3 ± 0.2. Then, the specific frequency results 
Sn = 1.3 ±0.15. 


3.5 NGC 4753 surface photometry 

NGC4753 presents several dust filaments crossing its centre, 
mainly in the East-West direction, as well as faint plums in the 
outskirts of the galaxy. Since their presence could point to a dis¬ 
torted structure, we obtained NGC 4753 surface brightness profiles 
in the Gemini GMOS-S g', r' and i' exposures and the MOSAIC II 
R image applying the task ELLIPSE, within Iraf. 

Figurel^shows the radial dependence of some important mor¬ 
phological parameters, derived for g', r' and i' filters. No striking 
difference is seen among them in ellipticity, position angle or the 
higher harmonic A4 associated to diskyness/boxyness. The strong 
isophote twisting which is present in the inner part of the galaxy 
(at radii smaller than 40 arcsec) is not shown since ellipses are not 
representative here. It is apparent that the variation in the PA at 
small radii is wider in g' than i, as expected from dusty substruc¬ 


Figure 9. From top to bottom, ellipticity, position angle, harmonic A4 and 
{g' — i')o colour profile in the radial regime 30 — 160 arcsec, obtained 
from our GMOS-S observations. In the three upper panels, green circles, red 
squares and black triangles represent parameters from g', r' and i' filters, 
respectively. 


tures. The parameter A4 is mainly positive, which is expected in 
disky isophotes, from rotationally supported galaxies. The {g'—i')o 
colour profile presents a slight gradient, becoming bluer towards 
the outskirts. 

The filled circles in the lower panel of Figure[To]represent the 
profile in the R filter. For comparison, the surface brightness pro- 
file in the r' filter w as overimposed, assuming r' — R = 0.25 
jFukugita et al.lfl99^ . Both profiles match within 200 arcsec. The 
differences for larger radii are due to the smaller GMOS field of 
view. The upper panel from FigurefTOl shows the residual map ob¬ 
tained from the subtraction of the R galaxy model to the origi¬ 
nal image. The image spans 12 x 12 arcmin^. There is clear ev¬ 
idence of underlying substructure. The most striking features are 
two opposite lobes, that seem to be connected to the galaxy centre 
by spiral remmants. Sparser spiral remmants are identified in the 
outskirts. 


4 DISCUSSION 

4.1 The paucity of GCs in galaxies with Type la Supernovae 


The work on Supernovae rates, specially over the last years has led 
to a current knowledge of this subject where sta r formation rates 
(hereafter SFR) are unarguably linked to SN rates l Cappellaro et al.l 
I 1999 I : iMannucci et alj2005l : ISullivan & et al.ll200^ . 

While it is natural to understand the presence of core-collapse 
SNe in starburst galaxies, the case of type la SNe in early-type 
galaxies is particularly intriguing. 
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Figure 10. Upper panel: Residuals map, obtained from the subtracion of 
the ELLIPSE galaxy model to the MOSAIC II R image. The image spans 
12 X 12 arcmin^. North is up, East to the left. Lower panel: Smooth lumi¬ 
nosity profile of NGC 4753 obtained from the MOSAIC II R (filled circles), 
and from the GMOS-S r' (open circles) images . The latter one was scaled 
assuming r' — R = 0.25 iFukugita et al.ll99^ . 


In their study of the correlation betwee n age, metallicity and 
SNIa peak magnitude. lGallagher et al.l ( l2008h find that the presence 
of SNIa in an elliptical galaxy is not dependent on its age or metal¬ 
licity. They claim, however, that peak la luminosities are correlated 
with age, while the trend observed with metallicity is an artifact. 
In addition, they observe that the rate of la SNe is higher in E/SO 
galaxies with spectroscopic ages younger than 3 Gyr than in their 
older counterparts. 

Given the mounting evidence favouring the SFR-SNela con¬ 
nection, one would expect to see young or intermediate-age pop¬ 
ulations in E/SO galaxies that hosted SNIa. Indeed, some striking 
examples exist of ear ly-type gal axies hosting SNIa’s (NGC 1316 
Stritzingeret^l2010l. NGC 4374 lCantiello et al.ll201 iL NGC 5128 
Jha et al j2007l) . 


In NGC 1316, interme diate-age globular c lusters have been 
spectroscopically suspected i Gomez£^l.l l200ih and spectroscop¬ 
ically confirmed jGoudfrooii et al.ll200ll) . Also the kinematics of 
this GCS is particularly complex, probably as a consequence of 
the 2 Gyr merger event jRichtler et al.ll20i4) . Moreover, the GCs 
colour distribution presents a distinctive composition, with a bright 
population in the bu lge that might be do minated by young and in¬ 
termediate age GCs jRichtler et al .l2012h . 


A common feature of these GCSs, which is observed in 
NGC 4753 as wel l , seems to be a low specific frequency (e.g. 
iGomez et alj|200i] : iRichtler et alj|201^ . This could be explained 
by a dearth of globulars or by a luminosity excess caused, for in¬ 
stance, by the presence of younger stellar populations. 


4.2 NGC 4753 and the general scheme of lenticulars 
evointion 


In the current scheme, the evolution from spirals to lenticular galax¬ 
ies in dense environments is due t o their interaction with the me¬ 
dia or with their close neighbours jMooreet^lJ 199^Bekkiet_S 
20021: ISun et al.l l2006l: iBarr et al.l l200l IWilman & ErwinI 12'^ 


Johnston et al.l2Ql'4 ). However, this mechanisms fail to explain the 


existence of SOs in low density environments, like the one where 
NGC 47 53 is located in, a n d alternative or igins should be con¬ 
sidered dStocke et al.ll2004l) . iTal et al.l (l2009h studied a sample of 
bright field Es and considered that the tidal features found in the 
diffuse light distribution of most of them, are evidence of mostly 
recent mergers. In that sense, the tidal features in NGC 4753 could 
also be related to a past minor merger. 

On the basis of a “mock” catalo gue created from the Millen¬ 
nium Simulation, iNiemi et"^ ( 1201 Oh have shown that almost all 
isolated E galaxies have undergone merging activity during their 
formation history. This is supported by the evidence of merger 
events in SO galaxies in low densi ty environments , e.g. t he case of 
NGC 6861 which was identified bv lEscudero et al.l ( l2015h by study¬ 
ing their GCs. Inhomogeneities in the projected spatial distribution 
of GCs, as spiral-like features, streams and/or arcs, have also been 
found in othe r early-type galaxies and inter preted as remnants of 
past mergers jP’Atosco et al.|[2013Ll2014al l3). On the other side, 
several galaxies that are likely remnants of recent mergers have 
GCs of intermediate-age (about 1 — 3 Gyr ), as shown by mean s 
of combined optical and NIR photometry by lTrancho et alj(l2014l) . 
That is, in different ways, GCs may be part of such intermediate- 

age populations. _ 

_ ^In the framework of the Project ATLAS®® dCappellari et al.l 

ItoTH) , a study of the distributi on and kinematics of CO gas in a 
sample of early-type galaxies jAlatalo et alj|201^ has been pre¬ 
sented, being NGC 4753 one of their targets. According to the 
CO maps, this galaxy is assigned a morphological classification 
M (“Mildly disrupted”) due to the intrincated dust filaments. Al¬ 
though generally dust and CO have similar distributions, NGC 4753 
is one of the clear cases in which dust is more extended than the 
detected CO (their Fig. 3). The CO velocity map of NGC 4753 is 
presented in their Fig. 2 and A20, showing that NGC 4753 is a 
massive regular rotator. On th e basis of the firm statement f rom 
ISteiman-Cameron et alj ( Il992h . regarding that the dust distribution 
may be understood as a processing disc, the molecular CO disc can 
be also considered unsettled, with an external origin (e.g. an accre¬ 
tion event). 

Unusual behaviours in the colour distribution (e.g . 


G6mez_et_alJ_ 2001; IRichtler et al.l 201^^ Caso et aP 2013bl : 
Escudero et al. l2015h . or GCLF (e.g. iGoudfrooii et al.ll2007l) ha\^ 


been considered evidence for GC subpopulations with intermedi¬ 
ate ages, asociated with late mergers. In some cases, it has been 
spect roscopically confirmed jGoudfrooii et al.|[200ll : leader et al.l 
l2004h . These features are also present in NGC 4753 GCS, pointing 
to the presence of GCs with intermediate ages, probably related to 
the same stellar population that produced the SNIa. In this sense, 
another hint in NGC 4753 is the spatial projected distribution of 
GCs, which differs from the usual scenario in early-type galaxies. 
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with red GCs concentrated tow ards the galaxy, whi l e blue ones 
form a more extended ha lo (e.e. lBassino et al.]|2006al : rUsher et alj 
l2013l:lDurrell et al.ll2014h . 


5 SUMMARY AND CONCLUSIONS 

In this work, the first detailed optical study of the GCS surround¬ 
ing the galaxy NGC 4753 is presented. This bright galaxy (Mv ~ 
—22.3), classified as lenticular as well as irregular, presents a re¬ 
markable dust structure with large filaments and is located in a 
low-density environmet. All these characteristics make it a good 
target to study galaxy evolution in the context of the current galaxy 
formation theories. 

On the basis of GMOS g', r', i' images obtained at Gemini- 
South Observatory, and Mosaic II R images from the Blanco 4-m 
telescope at Cerro Tololo Interamerican Observatory, we performed 
the photometry of the GC candidates and surface photometry of the 
galaxy itself. As expected for this galaxy, several properties of the 
GCS are peculiar. The GC colour distribution can be well fitted by 
a trimodal set of blue/intermediate/red subpopulations, instead of 
just the common blue/red one. The fraction of blue GC candidates 
is similar to that of intermediate plus red ones. The intermediate- 
colour subpopulation can be explained as a SSP of an intermediate 
age, i.e. 1 — 3 Gyr. The GCLF can be fitted by a Gaussian, but an 
excess of bright globulars is present for i' > 23 mag. If the pro¬ 
jected distribution of such a bright subsample is analyzed, it shows 
inhomogeneities in their azimuthal distribution, with fewer candi¬ 
dates towards the North. Integrating the GCLF, an estimation of 
1070 ± 120 GC candidates is obtained for the whole GCS, which 
leads to a specific frequency of Sn = 1.3 ± 0.15, also a quite low 
value for such a massive galaxy. 

The surface photometry of the galaxy reveals a complex un¬ 
derlying substructure. In addition to the twisted dust lanes and fila¬ 
ments, a pair of inner distorted spiral arms linked to another pair of 
outer ones appear when a smoothed light model is subtracted from 
the R image. Some plums in the outskirts are also visible in the 
residual map. Another peculiarity of NGC 4753 is that it hosted two 
SNIa (SN1965I and SN1983G), i.e. a high SNIa rate for a lenticular 
galaxy, that is usually related to the existence of intermediate-age 
population. We suggest that the intermediate-colour subpopulation 
of GCs is probably one of the components of that population. 

All these properties of NGC 4753 and of its GCS, converge to 
a scenario of a quite recent accretion event or merger, that might 
have taken place 1 — 3 Gyr ago. The intermediate-age population 
that might have formed then, may be identified by the intermediate 
colour GC subpopulation and the high SN rate of this galaxy. 
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